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Abstract. The PHENIX Run4 data-set provides a powerful opportunity for 
exploring the angular anisotropy of identified particle yields at high px- Com- 
plementing traditional v-i measurements, we present tt° yields as a function 
of emission angle with respect to the reaction plane in Au+Au collisions at 
^snn = 200 GcV/c. The centrality dependence of the angular anisotropy 
allows us to probe the density and path-length dependence of the energy loss 
of hard-scattered partons. 
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1. Introduction 

Two of the greatest mysteries that have arisen from the RHIC physics program are 
the source of the apparent flatness of the high pt (> 5 GeV/c) suppression of Raa [ 
^ and the source of non-zero V2 at high pt [H] ■ The existence of intermediate to high 
Pt V2 was suggested early in the RHIC program [E], and has been the subject of 
many theoretical treatments (see [0] E] for some additional examples). Traditional 
flow and parton energy loss pictures have failed to describe the magnitude of this 
anisotropy. Measurement of the azimuthal asymmetry t>2 at high px will shed 
light on the contributions from flow, recombination, and energy loss, as well as the 
transition from soft to hard production mechanisms. 



2. Measuring t> 2 and n yields in PHENIX 

The orientation of the reaction plane is measured event-by-event using the set of 
two PHENIX Beam-Beam Counters (BBCs), which reside at the region 3 < < 4. 
Each detector is an array of 64 hexagonal, close-packed quartz Cherenkov counters, 
located 150 cm from the interaction point. The charge measured by each counter 
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is proportional (on average) to the multiplicity of particles hitting it. The reaction 
plane angle ^pp is determined from the value of (cos2<j)). Because the two BBCs 
provide independent measurements of ^pp, we can estimate the resolution of the 
combined measurement via standard techniques 

For measuring photons and ir°s, we use the Electomagnetic Calorimeter (EM- 
Cal) [Cj. Candidate clusters are required to pass 7 identification cuts, and n%i nv 
distributions are formed from pairs of these clusters. The resulting yields arc binned 
in angle with respect to the reaction plane (Acj) = cj) — ^pp). A similarly binned 
mixed event background is then subtracted. The counts in the remaining peak cen- 
tered on the 7T mass are integrated in a ±2er window (where a is the width of a 
Gaussian fit to the peak). Six bins in Acj) are used in the interval [0 — 7r/2]. 

To measure v 2 , we fit the Acj) distribution as 1 + 2v 2 aw cos(2A0). The resulting 
v 2 parameter needs to be corrected for the reaction plane measurement resolution, 
hence the designation v™ w . The resolution opp is determined for each centrality 
bin, and leads to the corrected value v 2 orr = v 2 aw /app. The yields as a function of 
Acf> can then be corrected with a factor 

l + 2^r r co S 2A0 
f{A4)) = 1 + 2,— cos 2 A*" (1) 

3. Results and Discussion 

To obtain Raa(A4>), we exploit the fact that the ratio of the yield at a given Acfs 
to the inclusive yield is equivalent to the ratio of the angle-dependent Raa to the 
inclusive Raa- Thus multiplying these relative yields by an inclusive measured 
Raa, we have: 

Raa{&4>) = Yield(A</>) /Yield x R A a (2) 



The Raa(A(/),pt) as a function of {N part } is shown in Figure 1(a) It is clear 
that there is non-trivial angular substructure of the Raa, and that it varies with 
centrality. This feature is emphasized by plotting the data on a semi-log scale, 
showing that the Raa behaves differently in different A(f> bins. 



The resulting ir° v 2 is shown in Figure 1(b) For the first time we observe 
v 2 up to 10 GeV/c. While the value of the v 2 decreases beyond intermediate px, 
it nonetheless shows a substantial and perhaps constant value out to the highest 
measured transverse momenta. 

To gain insight into the v 2 mechanisms at work at high px, we turn to models. 
We compare the ir° v 2 to two models, a calculation done by Turbide et al. (using 
an Arnold-Moore-Yaffe (AMY) formalism [ 0) and the Molnar Parton Cascade 
(MPC) model [ HUj . Figure[21shows calculations from these models, plotted alongside 
data for similar centralities. The AMY calculation contains energy loss mechanisms 
only, and we see that the data appear to decrease to a value at high p T that is 
consistent with this model; the level of agreement is most striking in the 20-40% 
bin. 
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The MPC model has a number of mechanisms, including corona effects, energy 
loss, and the ability to boost lower px partons to higher px (a unique feature). The 
calculation shown in Figure [21 does a better job of reproducing the overall shape of 
the V2, though it is systematically low. It is important to note that this calculation 
is done for one set of parameters, so it should be very interesting to see if the MPC 
can better reproduce the data for a different set of parameters (the opacity of the 
medium, for example). 

The prevailing thought is that the high px behavior of the i>2 is due to energy 
loss mechanisms. If this is true, the Raa should be sensitive to the geometry of 
the collision. To test this behavior, we seek to reparamaterize the two handles we 
have on geometry (centrality, or collision overlap, and angle of emission) into a 
single parameter, a quantity which we will refer to as " p L dL" . Details of the 
calculation are described in [E]- This effective path length is calculated from the 
parton-density weighted average of the length from hard-scattering origin to edge 
of an ellipse and includes the Bjorkcn 1-D longitudinal expansion. We perform a 
Glauber Monte Carlo sampling of starting points to account for fluctuations in the 
location of the hard-scattering origin of the particles' paths within the region of 
overlap between the colliding nuclei. The crucial feature of p L dL is that it is 
proportional to the energy loss sustained by the parton as it traverses the medium. 

The result of plotting Raa for all centralities and angles vs. p L dL is shown in 
Figure |3 If the observed Raa arose from only geometric effects, we would expect 
the data to exhibit a universal dependence on p L dL. For low px, this is clearly 
not the case; something more than just energy loss is taking place there. However, 
when the px reaches 7 GeV/c and above, the Raa data do indeed appear to have a 
dependence on a single p L dL curve. This apparent scaling strongly suggests that 
the dominant effect on Raa at high-p-r is radiative energy loss. 

4. Conclusions 

We have presented the first measurement of high px «2 for ir°s. It is now clear 
that the V2 at high px does decrease but to a non-zero value. Comparison of v<z 
with models suggest that the dominant mechanism at work at high px is energy 
loss. In addition, we have presented the first measurement of ir° Raa as a function 
of angle with respect to the reaction plane. The Raa{^-<I>tPt) exhibits interesting 
angular substructure. Furthermore, when the Raa data are plotted as a function 
of an effective path length through the medium, the scaling that arises at high px 
also argues for energy loss as the dominant mechanism at work. 
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(b) TT° V2 vs. Pt 

Fig. 1. 7T° Raa and v 2 - |(a)| shows Raa(pt) as a function of (N part ) for each A</> 



bin; the panels are for different pj> bins, (b) Shows the ir° V2{pt)i with each panel 
corresponding to a centrality bin. The grey bands indicate the systematic error due 
to the reaction plane resolution correction. 
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Fig. 2. Comparison of 7r° i>2 with models. The top three panels show the AMY 
calculation with data for three centralities [IB]- The bottom two panels compare 
two centralities with the b = 8 fm calculation of the MPC model [I1U|. 
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Fig. 3. Raa(A(1),Pt) vs. p L dL. The panels correspond to different px ranges. 
The solid circles are the most peripheral events, while the solid stars are the most 
central events. 



